Susceptibilities and spin gaps of weakly-coupled spin ladders 
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We calculate the uniform and staggered susceptibilities of two-chain spin-^ Heisenberg ladders 
using Monte-Carlo simulations. We show that the gap extracted from the uniform susceptibility 
and the saturation value of the staggered susceptibility are independent of the sign of the inter- 
chain couphng J± in the asymptotic limit \ J±\/J —> 0. Furthermore, we examine the existence of 
logarithmic corrections to the linear scaling of the gap with \J±\. 
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Spin ladders are arrays of coupled spin chains and ex- 
hibit structures that interpolate between a single one- 
dimensional chain and the two-dimensional square lat- 
tice. Their properties are peculiar as the magnetic spec- 
trum is gapless only in the case of non-integer-spin lad- 
ders formed of an even number of chains Spin-i lad- 
ders are thus closely related to spin-S chains for which 
the excitations are gapped for integer spins and gapless 
otherwise.'' An important question is the understand- 
ing of the transition between the gapless chain and the 
gapped ladder (the Haldanc phase) . That issue has been 
studied both analytically'*-^ and numerically!^ 

In this paper, we report results of Monte-Carlo simula- 
tions on two-chain spin-^ ladders with antiferromagnetic 
coupling along the chain direction and ferromagnetic or 
antiferromagnetic coupling across the rungs. We calcu- 
late the uniform and staggered susceptibilities and show 
that, in the asymptotically weak rung-coupling regime, 
these thermodynamic quantities are independent of the 
sign of the rung coupling. In this regime, our results 
for the spin gap are consistent with the theoretically- 
predicted existence of logarithmic corrections to the lin- 
ear behavior A ~ | ^± | • 

The spin-i Heisenberg ladder magnet is described by 
the following Hamiltonian: 



N 



N 



i-\-l,n 



(1) 



= 1 n=l 



where i runs along the chains, N is the length of the 
chains, J is the antiferromagnetic coupling along the 
chain direction and J± is the coupling between the two 
chains (the rungs of the ladder). We consider both anti- 
ferromagnetic {J± > 0) and ferromagnetic {J± < 0) rung 
couplings. The Hamiltonian Eq. (1) is investigated with 
the loop cluster algorithm with a discrete Euclidean time 
grid.®^ Periodic boundary conditions are used along the 
chain direction as well as in the Euclidean time direc- 
tion. The chain length is kept larger than six times the 
spin-spin correlation length^. The Trotter number, which 
defines the discretization of the Euclidean time axis, is 
kept larger or equal to 20 /T where T is the temperature. 
We calculate the uniform and staggered susceptibilities. 
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FIG. 1: Uniform (top) and staggered (bottom) susceptibili- 
ties as function of the inverse temperature for the spin-^ lad- 
der. Left (right) panels are for rung couplings \ J±\ — 0.10 J 
(|Jx| = O.OIJ). 



defined as: 
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Studies at intermediate and large couplings show that 
finite-size effects, both in the Euclidean time direction 
and in the chain direction, are smaller than the statisti- 
cal uncertainties for the two observables when the lattice 
dimensions were kept larger than the previously defined 
limits. The same relations for the required system size 
were assumed to apply in the small coupling regime. 

Figure ^ shows examples of the uniform susceptibility 
and staggered susceptibility for ferromagnetic and anti- 
ferromagnetic J_L. For \J±\ = O.OIJ we find that the 
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FIG. 2: Low-temperature staggered susceptibility as a func- 
tion of the rung coupling. 
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FIG. 3: Spin gap as function of the rung coupling for both 
ferromagnetic and antiferromagnetic rung couplings. The 
dashed lines are fits to D [ln{\J±\/ J)]^ with exponent E of 
0.5 and 0.23 for, respectively, for ferromagnetic and antiferro- 
magnetic rung coupling (Eq. ® and Eq. JTJ, respectively). 
The continuous lines are fits to Eq. 0, with coefficients 
A = 0.53(2), B = 2.2(5), as well as C = 3.9(5) (C = 2.0(6)) 
in the antiferromagnetic (ferromagnetic) case. 



susceptibilities become independent of the sign of J±, 
consistent with theoretical expectations for two weakly- 
coupled Heisenberg chains.^ 

Figure |21 shows the dependence of the extrapolated 
zero-temperature staggered susceptibility on J±/J. In 
the weak-rung-coupling regime, Xs(0) varies linearly with 
\J/J±\, as can be seen in the bottom panel of Fig. [3 We 
find that \ J±\xsiO) = 11.5(5) for \ J±\/J 0. 

The spin gap is extracted from the uniform suscepti- 
bility using the following low temperature form:** 



-l/2p-A(Ji)/T 
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where A is the gap. This limit is derived from a non- 
interacting magnon model with a quadratic band dis- 
persion. The fits are performed for T < A/3, and the 
results are shown in Fig. |31 The data for antiferromag- 
netic rung couplings { J± > 0) agree well with previously 
reported results.^ In the weak coupling limit, the gap 
deviates from a linear dependence on J±. For the ferro- 
magnetic rung coupling, values of the gap at intermediate 
rung couplings (0.1 < \ Jj_\/J < 1) are in agreement with 
previous reportsiSiii*i^. 

At very large ferromagnetic rung coupling, the spin- 
^ ladder is equivalent to a spin-1 chain with coupling 
J/2. The gap for the spin-1 chain was calculated to be 
A = 0.41050(2)^ or A = 0.4107(1)J" by the den- 
sity matrix renormalization group (DMRG) technique. 
In the present investigation, the spin-1 chain gap was es- 
timated at 0.406(3) J, which agrees reasonably well with 
the DMRG value. At finite ferromagnetic coupling J±, 
the ratio A/J± increases with decreasing value of J±. 
At small rung coupling, the data for the antiferromag- 
netic and ferromagnetic rung coupling exhibit a similar 
trend. Together with the behavior of the staggered sus- 
ceptibility, this seems to confirm the assertionSiii that 
the properties of a two-chain spin- 1/2 spin ladder are 
independent of the sign of J± in the weak coupling limit. 

In the weak-coupling regime, the gap has a nearly lin- 
ear dependence on the rung coupling J±. Logarithmic 
corrections have been proposed to this linear dependence 
based on field theoretical considerations. Shelton, Ners- 
esyan and Tsvelik^ concluded that the gap follows the 
form: 
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where A, B and C are unknown constants. Totsuka and 
SuzukiiS, on the other hand, suggested that the gap be 
described by 
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For |J±|/J < 0.1, Eq. ^ describes the gap well 
for both antiferromagnetic and ferromagnetic rung cou- 
plings, while Eq. © satisfactorily describes only the 
ferromagnetic result. A modified form of Eq. © is nec- 
essary to obtain an acceptable description for antiferro- 
magnetic rung couplings: 
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Although Eq. (5) describes both results well below 
\Ji_\ = 10" V, it would be necessary to know the gap at 
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rung couplings well below |Jj^| — 10^^ J in order to con- 
clusively establish that this form is asymptotically correct 
in the limit | Jj_|/J — > 0. 

In conclusion, we have calculated the uniform and stag- 
gered susceptibilities for the spin-i two-chain Heisenberg 
ladder. We establish that the susceptibilities are inde- 
pendent of the sign of Jl in the weak-coupling regime 
\J\-\I J ^ 1- III that regime, the staggered susceptibility 
is linear in J/| Jj^|, while the excitation gap is described 



by a logarithmic correction to the linear dependence on 
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